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SUMMARY 


During an eight-week visit approximately 200 simulations using the 
computer model for the regulation of ery thropolesls were carried out In four 
general areas: with the human model simulating hypoxia and dehydration, 

evaluation of and simulation of dehydration using the mouse model. The ex- 
periments led to two new considerations for the model. Firstly, a direct 
relationship between erythropoietin concentration and bone marrow sensitivity 
to the hormone and, secondly, a partial correction of tissue hypoxia prior to 
compensation by an increased hematocrit. This latter change in particular 
produced a better simulation of the effects of hypoxia on plasma erythropoietin 
concentrations. The model produced excellent qualitative changes of dehydra- 
tion. In addition, a greater appreciation of the gaps in our knowledge of 
the regulation of erythropoiesis became evident. There was a considerable 
cross-fertilization of ideas and concepts between the author and the computer 
group. As a result, at least four publications are in the formulation stage. 


INTRODUCTION 


A computer model for the study of erythropoietic control mechanisms 
with particular emphasis on the "anemia" of space flight has been developed 
and described by Leonard (1,2,3). The model is capable of simulating the 
effects of various erythropoietic stimuli on red cell production. However, 
a full appreciation of the model was difficult as it was developed by personnel 
with principally a mathematical background with little input from biologists. 

The primary purpose of the present author's leave of absence from the University 
of Tennessee to The Johnson Space Center was to interact with the computer 
groups with the aim of further validating the model of erythropoietic control 
mechanisms and to determine how well a modified model for mouse parameters 
responded to data obtained from a potential animal model of the "anemia" of 
space flight. As a secondary objective the author was asked to advise in the 
most appropriate role, from a grantee's viewpoint, for NASA scientists in the 
Spacelab Life-Sciences Program. The leave of absence extended from 14 August 
to 5 October 1978. 

THE MODEL 

A detailed description of the General Electric model for the control 
of erythropoiesis can be found in the reports by Leonard (1,2,3) and Grounds (4). 
Salient points only will he reviewed in this report. The model is summarized 
in Figure 1 (from reference 1). It consists of 69 parameters of which 24 are 
variable, the remainder being changed automatically in response to the input. 

The model was established using numerical values appropriate for the human 
and it has been shown to realistically simulate a variety of erythroid 
perturbations eg: compensated hemolytic anemia, r?d cell infusion (2). Hypoxia 

and the "anemia" of space flight cannot be directly simulated by the model (!•). 












RESULTS AND DISCUSSION 


Approximately 200 simulations were carried out. These can be grouped 
into four categories: 

1. Simulation of hypoxia using the human model, 

2. Simulation of dehydration using the human model but with data 
"scaled-up" from studies In mice, 

3. Evaluation of a model established with numerical values 
appropriate for mice, 

4. Simulation of dehydration using the mouse model. 

I. Sim ulation of hypoxia using the hu man model: 

The model responds to 7 days of hypoxia with an acceptable change in 
hematocrit (HCT, Figure 2) and red cell volume (RCV, data not shown). However 
it predicts that erythropoietin production (EP) and plasma erythropoietin con- 
centration (ERY) reach peak levels by the third day and are maintained at these 
high levels for the duration of hypoxia (Figure 2). In contrast, experimental 
data suggest that plasma erythropoietin concentrations respond to hypoxia by 
reaching a peak at 1-3 days and then fall abruptly to titers only slightly 
above normal despite continuation of the stress (5-9). The object of the 
initial simulations was to determine the affect on the HCT of allowing ERY to 
reach maximum values and then reducing it to a level approximately twice 
control. ERY has altered by changing EPO (Mnemonic II 49); the extrapolated 
value of EP at zero tissue oxygenation. Simulations determined that if EPO 
was altered prior to day 3 of hypoxia the maximal level of ERY was not reached 
and that a reduction of EPO from 20.09 to 8.0 was necessary to reduce ERY to 
approximately twice normal. The effect of changing EPO at day 3 in a 7 day 
exposure to hypoxia is shown in Figure 3. As might have been anticipated, such 
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FIGURE 3 
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a change resulted in a smaller increment In HCT than was observed when ERY 
wan allowed to remain high (compare Figures 2 and 3). This difference became 
more marked if the duration of hypoxia was extended to 21 days when the HCT 
only reached 46.87! (Figure 5, compared to 50. 9X when ERY was allowed to remain 
high (Figure 4). Further simulations were done in an attempt to bring the HCT 
up to levels observed when ERY was allowed to remain high. The rationale for 
the following simulations was based on recent studies using in vitro clonal 
assays for Erythropoietin Responsive Cells (ERC). These investigations (10,11), 
together with earlier computer simulations (12), have suggested a direct 
relationship between ERY and the numbet and/or sensitivity of the ERC population. 
TMs change in number and/or sensitivity is best expressed in the present model 
by a change in G2; viz, the Cain factor relating erythropoietin production to 
red cell production. Adjusting G2 from 2 to 4 at the end of day 1 of hypoxia 
and bringing ERY down to twice normal on day 3 (Figure 6) results in increases 
in HCT, RCV and Blood Volume (BV) comparable to that seen (Figure 4) when ERY 
is allowed to remain at values approximately six times normal for the duration 
of exposure to the hypoxia stress. Other simulations (not shown) determined 
that a doubling of G2 on day 1 was the minimum change necessary to elecit the 
necessary increment in HCT although a similar response could be obtained by 
altering PI - a comparable parameter to EPO but relating to the bone marrow 
response (Figure 7). Subsequent studies (not shown) demonstrated that G2 
could change in parallel with ERY or could remain high throughout hypoxia with 
little effect on the HCT at the end of 21 days of hypoxia; G2 could, therefore, 
be directly related to ERY. Partial confirmation of thi3 hypothesis was 
obtained by discussions with Dr. C. Peschle at the annual meeting of the 
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IMULATION OF HYPOXIA (P02 ART = 55.) WITH CHANGES IN 62 (DAY 1) AND EPO (DAY 3) 
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Simulation of hypoxia (P02 ART = 55.) with changes in PI (day 1) and EPO (day 5). 
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International Society for Experimental Hematology. His, currently unpublished, 
data suggest that there are tempoiary net Increases in ERC In response to 
hypoxia although these cells In the bone marrow and spleen behave somewhat 
differently. There Is, therefore, some evidence that ERY can directly influence 
a factor that may be mathematically designated G2. However, at this time the 
biological mechanism whereby ERY would decrease In the face of maintained 
hypoxia was uncertain. 

In order to more fullv understand the relationship between G2 and FRY 
changes post-hypoxia were simulated with (Figure 8) or without (Figure 9) changes 
in G2 and a reduction in EPO. There was little evidence that changes in G2 
and EPO had a detrimental effect on the simulation of the post-hypov-'a changes; 
in fact, it was in the situation where G2 and EPO were altered (Figure 8) that 
a more realistic post-hypoxia suppression of P was obtained. 

Ir* conclusion to this part of the study it appeared that if G2 was 
altered in a way which reflected changes in ERY, then ERY could be brought 
back to titers aporoximatcly twice control without a detrimental effect on 
red cell production. The original i.iodel did not consider this possibility or 
a mechanisr .or bringing FRY down from its peak levels. Two modifications 
were made to the .nodel : 

ti) a direct i 'urical relationship between ERY and G2 was 
incorporated Into the Subroutine Blood, i.e., 

G2 - ERY x 2 

(ii) Since the only 'trigger' in the model for changes in erythropoietin 
production is tissue oxygenation, it was presumed that a mechanism 

• « « « • • i 4 

• • 4 • 4 « « 4 

would be operative whereby tissue oxygenation would be brought 
back almost to control levels Before the i matocrit completely 
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FIGURE 9 



Simulation of m days of hypoxia (P02 ART = 55.) and post-hypoxic changes in the 

HUMAN MODEL. 
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compensated for the hypoxia. It was postulated that such a 

h 

mechanism could occur by changes in blood flow, rate of oxygen 
diffusion, or in the oxygen-carrying capacity of the blood i.e., 
the p50. On this basis the following formula was inserted into 
Subroutine Tissue, line 24: 

RDIF02 - RDIF02 + (P02TIS } /SENS - RDIF02)/DAMPF 
Where RDIF02 ■ rate of oxygen diffusion, RDIF02 ■ original rate 
of oxygen diffusion, P02TIS ■ level of tissue oxygenation and 
SENS and DAMPK ■ numerical values to alter the level of sensitivity 
of the new model. In addition a delay parameter, DLAY, was 
introduced to allow a variable period of time to occur before 
partial correction of the tissue hypoxia was instituted. 

These modifications were incorporated into a modified model 
designated "0 ADD THERM. LOADCD." Nineteen simulations were per- 
formed to optimize the 3 new variables at 

DLAY ■ 2880 minutes 
DAMPF ■ 5760 minutes 
SENS - 0.55 

and the simulation of hypoxia shown in Figure 10 was obtained with 
this model. Note the changes in HCT, RCF and BV compere well with 
the original model (Figure 4) in which ERY was allowed to remain 
high. Note also the partial correction of P02TIS which accounts 
for the decrease in ERY (and, therefore, G2) . Preliminary experi- 
ments in this laboratory using mice exposed to hypoxia suggest 
that there are dramatic shifts in the p50 which might biologically 
explain the partial correction of P02TIS. "Dose/reeponse" 


FIGURE 10 
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relationships to various degrees of hypoxia and the changes during 
the following 14 days post-hypoxla are summarized In Figures 11 
through 14. At hypoxia levels of P02ART t 35 the system aborts. 

2. Simulation of dehydration using t he human model but with data "scaled-up" 
from studie s In mice : 

We (13,14) have obtained much Information regarding the effects of 
dehydration and/ot food consumption on erythropolesls mice. Thus we have 
data on changes In plasma volume (PV), ERY, In the sensitivity of the target 
tissues to erythropoietin (G2) , in red cell production (P) and an Indication 
of changes In metabolic rate (AOM) assumed to be reflected by a change in body 
weight. By changing PV, G2 and AOM as die it*»d by experimental findings we 
hoped to simulate the changes In P and ERY. Since the data had been accumulated 
In mice, results expressed as percentage changes were used In the human model. 

(1) Simulations in the original model: 

The effects of programmed changes In PV, G2 and AOM on P and 
ERY over a three day period are shown in Figure 15. Note that in 
contrast to the experimental findings (13) ERY was decreasing at 
the end of day 1. In addition, suppression of P was much less 
than was observed experimentally (13,14). Suppression of P to 
values comparable to those seen In the mice at 24 hours could be 
accomplished by reducing the parameter (Z) for bone marrow cell 
transit time from 5760 minutes to 1 minute (Table 1). ERY can 
be maintained above 0.9 of control by reducing renal blood flow 
(BF) from 1.2 to 1.0 (Table 2) but under these circumstances P 
is higher than expected from the experimental /lata. . A three !oy 


0 « 


FIGURE 11 



MfcP D COPY WANTED' Y,N>. . . TIME ' DAYS 









FIGURE 12 







FIGURE 13 


dRt&NAL PA<fc * 
OF POOR QUALITY 



a 

o 

z: 


< 

x 

3 


a 

o 


LLI 

X 


X 

0 

CL 

> 

X 

1 

H- 

CO 

o 

a. 

a 

< 


LH 

LTl 

II 


OC 

<c 


CNI 

O 

Q_ 


X 

o 

a. 

>- 


3 

X 









I 


cfatGlNAL PA<jE Rk 
OF POOR QUALITY 


I 


! 




a 

o 

x 


z 

I 

3 

X 

a 

UJ 


a 

o 


< 

X 

o 

CL 

>- 


CO 

o 

o. 

a 

z 

< 


LT\ 

cr 

II 


cn 

<c 


CNI 

O 

Q_ 


< 

X 

o 

CL 

> 

X 


o 

z 

o 


3 

X 


f 

1 


V 


CO 








FIGURE 15 



Simulation of dehydration by changes in PV, G2 and AOK using the original human model 
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(ill) Simulation of rehydratlon with the original model: 

Following three days of experimentally-dictated changes in 
PV, G2 and AOM (Figure 16) these parameters were returned to 
control values (rehydration) and the model run for a further 
seven days. The simulation is listed in Table 4. An excellent 
simulation of the experimental data (14) was obtained even to 
the slight decrease in KCT and small Increase in P. 

3. Evalu a tlon of a model established with numerical values appropriate for mice: 
A mouse model was established with the designation "0 ADD THERM. I.OADM." 
This model responded appropriately to hypoxia (P02ART ■ 55. or 35.) but 
inappropriately to Infusions of erythropoietin (ENIF) , to changes in e ythro- 
poietin production (ERYO), erythropoietin half-life ir.HL), infusions of red cells 
(RCFLO), changes in BF, PV, red cell destruction (RKC), basal metabolic rate 
(BMOTIS), and hemoglobin oxygen-carrying capacity (CHLQ2). These inappropriate 
responses were eventually traced to the number of significant figures used in the 
integration steps. When these were increased the model responded appropriately 
to the above stimuli in 18 simulations. As an example of the behavior of this 
model. Figure 18 depicts simulation of a typical erythropoietin assay similar 
to those carried out routinely in our laboratory. Mice were presumed to be 
exposed continuously to bl oxygen (P02ART ■ 25) for 14 davs before return to 
an ambient environment (P02ART - 78). On the seventh day after return 
the animals wete infused with erythropoietin at a dose twenty times the 
normal, pre-hypoxia, plasma concentration and the model run for a further four 
days. A qualitatively good simulation of the changes which are thought to occur 
during hypoxia and post-hypoxia were obtained. However, changes in PV during 
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hypoxia may occur. Note the small Increases in ERY and P in response to ENIF 
in comparison to those produced by exposure to hypoxia. Note also that P does 
not fail as much during the post-hypoxic period as experimental studies with 
radioiron would lead us to expect. 

As part of the validation of this model the balance between the 
sensitivity of the erythropoietin producing mechanism (Cl) and that of the bone 
marrow response to erythropoietin (G2) was elucidated. It was found that 
numerical values of G1 - 2.5 and G2 * 2.4 provided the most satisfactory increases 
in ERY and P in response to hypoxia. 

4. Simulation of dehydration using the mouse model: 

As detailed in Section 2 (above) we utilized our knowledge of changes In 
PV, C2 and AOM to determine how the model, specifically ERY and P, would respond 
to dehydration. The mouse model responded to the results obtained with mice 
in a similar way to the human model (Figures 15 and 16) using proportional 
changes from animal data. Like in the human model, optimum simulation of de- 
hydration (Figure 19) was obtained only if the bone marrow maturation time (Z) 
and renal blood flow (BF) were both reduced to 1.0 and 0.0015 respectively. 

Even under these conditions, red cell production (P) and red cell mass (RCM) did 
not fall to the low levels seen experimentally. Similar to the trials with the 
human model the HCT was subnormal and P slightly increased during the excellent 
simulation of rehydration — during which time Z and BF were returned to control 
values (Table 5). 

The reasoning behind the implementation of changes in Z and BF to get 
acceptable simulations of dehydration are not clearly understood. A decrease 
in blood volume (BV) might simplistical ly be expected to decrease BF unless 
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other compenaat lng mechanisms eg. heart rate, were Increased. On the other 
hand, decreasing the time (Z) for a cell to mature through the bone marrow 
from 5040 minutes to 1.0 minute Is probably unrealistic. A possible explana- 
tion of this dichotomy might be that dehydration Inhibits cell proliferation 
at a stage late In the maturation sequence. A reduction of cell output from 
the marrow would, therefore, occur earlier than when the Inhibition occurred 
at a more primitive cell level. However, even this theory Is not entirely 
compatible with experimental data as it would suggest changes in the marrow 
morphology during dehydration. No such change has been observed (14). 

On the basis of the above simulations and discussions of a theoretical 
nature with Drs. Leonard and Klmzey, four papers are In the formulation stage. 

CONCLUSIONS AND RECOMMENDATIONS 

The visit to the NASA group was a thought-provoking experience. It has 
led to a considerable "cross-fertilization" of Ideas between people trained in 
different disciplines; i.e., between a mathematically-oriented group and the 
author whose major experience is in the biological sclencee. As an example 
of this greater appreciation of our common interests in the regulation of 
erythropoieses, 1 was able to provide to the computer group an idea of the 
accuracy of the numerous assays involved and an indication of whether computer 
simulated data might be significantly different from those generated biologi- 
cally. In addition, several areas in which the model might be improved were 
recognized; e.g., relating ERY to G2, rapidly correcting for tissue hypoxia by 
mechanisms which we are currently attempting to elucidate. In return, it 
became evident from the simulations performed and from discussions with the 
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computer group that an understanding of the regulation of erythropolesla la 
not as complete is we had all thought. This has led to a greater appreciation 
of erythropoietin control mechanisms and underlined the gaps in our knowledge, 
some of which could be filled by experiments in my laboratory. These ideas 
have been generated as a direct consequence of the visit to NASA. Certainly, 
a clearer understanding of the erythropoietic effects of hypoxia and dehydra- 
tion were gained and several avenues of research which could be performed in 
our potential animal model for the "anemia" of space flight (13,14) were 
recognized. These experiments will be enhanced by access to the computer 
model now established on The University of Tennessee's computer system — 
another direct "spin-off" from the visit to Houston. The planned Joint publi- 
cations using the model will also act to the mutual advantage of both groups. 

On the basis of my visit It Is recoranended that the concept of having 
an outside scientist to work with the computer group for a couple of months 
per year be continued. This will allow different inputs into the erythropoietic 
control model with the result that more and more realistic simulations are 
obtained. In addition, as more and more scientists use the model It will 
certainly become more widely known and probably appreciated to a greater degree. 
There seems little doubt that working in parallel between the computer model 
and the laboratory bench will greatly enhance our understanding of the regulation 
of erythropoiesis. 

Three improvements are suggested which may make subsequent visits by 
other scientists even more useful: (a) the visiting scientist should attempt 

to spend one or two days with the computer group perhaps a month before the 
planned leave of absence to familiarize himself with the computer systems; 
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(b) aC lease one member of the computer group be allowed dt ^dcZo leave from 
his usual responsibilities and given the sole job of Interacting with the visitor. 
These two suggestions apply particularly to visiting scientists whose knowledge 
of computer technology Is severely limited; (c) the possibility of "return 
visits" by one or more of the computer group to the visiting scientist's 
laboratory — obviously there are many factors which would need to be clarified 
before this could be made operative. 

By way of summing up the two months' stay at the NASA facilities in 
Houston, the visit was very useful and thought-provoking. It has led to a 
greater appreciation of the gaps In our knowledge regarding erythropoietic 
control mechanisms and identified areas of future research both for the computer 
group and author. I am exceedingly grateful to friends and colleagues both In 
Knoxville and Houston who made the visit possible in the first place and for 
making It so pleasant and successful. 
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